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Abstract 

We prove rigorously the occurrence of zero-mode Bose-Einstein condensation for a 
class of continuous homogeneous systems of boson particles with superstable inter- 
actions. This is the first example of a translation invariant continuous Bose-system, 
where the existence of the Bose-Einstein condensation is proved rigorously for the 
case of non-trivial two-body particle interactions, provided there is a large enough 
one-particle excitations spectral gap. The idea of proof consists of comparing the 
system with specially tuned soluble models. 
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1 Introduction 



Rigorous proofs of Bose-Einstein condensation (BEC) for realistic systems is an ongoing 
challenge for almost eight decades with renewed interest after the successful experiments 
with trapped gases of alkali metals. In particular establishing a proof of condensation for a 
system of interacting particles turns out to be a real hard problem. However it is believed 
that this the real issue, in view of its experimental realisation in superfluid 4 He. 

For this work, we were inspired by the recent results on condensation for trapped 
gases, i.e. inhomogeneous systems because of the external trapping fields [1] and by the 
result on Bose-Einstein condensation for systems with a gap in the one-particle excitation 
spectrum with a Van der Waals family of two-body potentials by Buffet, de Smedt and 
Pule [2]. Although these two exact results are a priori unrelated, the first one being for 
inhomogeneous systems and the second one for homogeneous systems, both start from 
systems with a gap in the one-particle excitations spectrum. Heuristically, in both cases 
one might start from a system with a gap in the spectrum and with condensation in the 
groundstate. Contrary to systems without a gap, one can anticipate that the condensation 
for systems with a gap is stable under small perturbations. One might understand that 
condensation in the groundstate, which is energetically isolated by a gap, can survive the 
switching-on of a gentle interaction, and that fluctuations must be of macroscopic size to 
overcome this gap and lift particles out of the isolated groundstate. 

In this note, we turn our attention to homogeneous continuous systems (i.e. without 
external trapping fields) of interacting bosons in the standard thermodynamic limit. We 
prove that there is Bose-Einstein condensation for high enough density and low enough 
temperature, provided that there is a large enough gap in the one-particle excitations 
spectrum. To our knowledge, this is the first example of a proof of condensation for 
systems with general two-body potentials. 

We remark that for the trapped systems one proves in [1] in fact macroscopic occupation 
of the ground state, which is the main property of BEC. However, as these systems are not 
homogeneous in any sense, there is no notion of phase transition present. The presence of a 
trap yields a fortiori a discrete one-particle spectrum. In our models we also assume a gap 
in the one-particle spectrum, and we prove standard BEC including the phase transition 
(gauge symmetry breaking, off diagonal long range order,. . . ) which follows directly from 
the homogeneity of the system. Therefore we consider our result as a step forward to the 
realisation of a bridge between the homogeneous and the trapped case, leading to a concept 
of phase transition also for the trapped case. In particular it might link the notion of the 
standard thermodynamic limit and the so-called Gross-Pitaevskii limit which is used in 
the trap case. 

Another challenging problem posed by our results is of course: can one close the gap ? 
i.e. can one prove that BEC persists if one takes the limit of the gap tending to zero. 

These were our starting ideas. First we define our systems under consideration. Let 
us consider a gas of interacting bosons in hypercubic boxes A = [-L/2, L/2] u C H u of 
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dimension v > 1 with periodic boundary conditions. The generalisation to other shapes 
for A is however readily done. Denote by V = L v the volume of the box A. The Hamiltonian 
for the volume A of this system in the boson Fock-space Tb reads 

Hl g = T£-fiN A + gU A , g>0, (1.1) 

where T£ is the kinetic energy operator with the gap A > in its spectrum, 

/I 2 A; 2 

t a =Z2 £ k a W - Aajao, e k = (1.2) 

We take units h 2 /2m = 1, the sum k runs over the set A*, dual to A, i.e. 

A* = {k G R u ; k a = 2irn a /L; n a = 0, ±1, . . . ; a = 1, 2, . . . , is} . 

The operators a k and a k are the Bose creation and annihilation operators for mode k G A*, 
the number operators are denoted by = a\ak, and the total number operator in the 
volume Ahy N A = EfceA* N k- 

We assume a priori the presence of a gap A in the one-particle excitations spectrum, 
isolating the lowest energy level. This might seem rather artificial, however note that the 
presence of such a gap can be realised in several ways. A gap can be created using attractive 
boundary conditions [3-5], i.e. considering systems in containers A with sticky boundary 
conditions. But we have to notice that the Bose condensate is not homogeneous in this 
case. Another possibility is to assume that part of the particle interaction has caused this 
gap and is as such effectively incorporated [6]. The idea of considering a gap is not new. 
In his book [7] F. London attempted to introduce the gap on heuristic grounds to clarify 
some of the spectral properties of superfluid 4 He. 

The interaction between the particles is modelled by the two-body interaction operator 
U\ = \ ( daxh/ a) (x)o) (y)v(x - y)a(y)a(x), (1.3) 

where oJ(x), o){y) and a(y),a(x) are the creation and annihilation operators for the Bose 
particles at x, y G H u . Below we assume that the pair interaction potential v(x) verifies 
the following conditions: 

(a) v : H u — > R 1 is a real, positive-type function from L 1 (R I/ ). As it is known (the 
Bochner theorem [8]), a continuous function is of positive type if and only if it is the 
Fourier transform of a positive measure d/i of finite total mass on H u : 

v(x) = / d/i(g) e iqx = v(-x). 
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(b) Since v G L 1 (R !/ ), the Fourier transform v(q) exists, and we suppose that 



v 



(0) = dx v(x) > and, v(0) > %), Vg G R", 



(1.4) 



where v(q) > by (a). 
Notice that by virtue of (a) the potential t>(x) at x = is finite and that v(0) > v(x), 



It is shown [9] that under conditions (a) and (b) the interaction is superstable, i.e. the 
n-body potential satisfies the inequality 



for some constants A > 0, B > 0, for all n G N, x% G A and A large enough which 
implies that the thermodynamic potentials exist for all values of the chemical potential 
fi. Moreover, our proof requires a more stringent condition: it asks the constant A = 
v(0)(l — e), with e > arbitrarily small and B = v(0)/2. These optimal stability constants 
for the L 1 -integrable potentials of positive type were established by Lewis, Pule and de 
Smedt [10]. Since these optimal constants are important for our proof, we incorporate in 
this note the corresponding argument (see section 2.1). 

The superstability of the particle interaction is, together with the assumption of a 
spectral gap (1.2), the physical foundation of our proof. We prove that the zero-mode 
Bose-Einstein condensation in the groundstate of the one-particle spectrum isolated from 
the continuous spectrum by a gap A is robust enough under the application of a superstable 
interaction, provided that this spectral gap is large enough. Therefore, we show that in 
this case the zero-mode BEC is stable with respect to a class of superstable interactions. 
More technically our main result, announced in [11], can be stated as follows: 

Theorem 1.1. Consider a system of interacting Bose particles (1.1) in three or more 
dimensions, with a two-body interaction satisfying conditions (a) and (b). Fix an inverse 
temperature and a chemical potential (/3,fj) such that fi > gv(0)p^((3), then there exists 
a minimal value for the gap A min , such that for all A > A min the thermodynamic limit 
(liniA : V — > oo) of the k = mode particle number occupation density is positive: 



(—)h£ (A/-0 denotes the grand- canonical Gibbs state for the model H^ g (1.1) at inverse 

temperature (3 = 1/ksT and chemical potential \x. v(0) is as in (1.4) and Pc(P) is the 
critical density of the perfect Bose gas at inverse temperature f3. 



Vx G W. 




Bn 



l<i<j<n 
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A detailed proof is found in section 3, Figures 1-2 illustrating this theorem can be 
found in section 4. Notice that an analogous result as in Theorem 1.1 holds in lower 
dimensions v < 2 (cf. section 4). Our proof of BEC is based on comparing the condensate 
density of the full model (1.1) with the condensate density of specially tuned reference 
systems. This yields (various) lower bounds on the condensate of the full system in terms 
of thermodynamic potentials of the reference systems. The reference systems we use are 
based on the mean-field Bose gas. Its thermodynamic properties, necessary in our proof of 
BEC, are reviewed in section 2.2. 

This paper is organised as follows: in section 2 the essential preliminaries, i.e. su- 
perstate potentials and the thermodynamics of the mean-field Bose g reviewed, 
section 3 contains the actual proofs of the new results and a discussion on these results 
(section 4) concludes the paper. 



2 Preliminaries 

2.1 Superstable Potentials 

An important aspect of particle interactions in both classical and quantum continuous 
systems are their stability properties [6,9,12]. If the particle interaction is such that the 
grand-canonical partition function does not converge, it is called catastrophic. It means 
that the thermodynamic pressure is not everywhere well defined, and good thermodynamic 
behaviour is excluded in those domains. 

A useful criterion for stability is the following: a pair-potential v is called stable [9] if 
the corresponding n-particle interactions can be estimated from below as : 

v(xi - Xj) > -Bn, (2.1) 

l<i<j<n 

for a B > 0, and all n > 2, Xi G A. This is a sufficient condition for good thermodynamic 
behaviour. An important subclass of the stable interactions are the so-called superstable 
pair-potentials [9], satisfying 

E v( Xi - Xj ) > ^n 2 - Bn (2.2) 

l<i<j<n 

for some constants A > 0, B > 0, and all n > 2, Xi G A, where A is large enough. This 
condition (2.2) yields the existence of the grand-canonical pressure for all values of the 
chemical potential ji G R. From (2.2), it follows that the interaction term (1.3) satisfies 
the operator inequalty : 

U A > 4j n I - bn a- (2-3) 
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An example of superstable potentials are those verifying the conditions (a) and (b). For 
this class of potentials Lewis, Pule, and de Smedt [10] proved the existence of optimal 
constants A = -0(0) (1 — e) and B = v(0)/2 in (2.2). Here, e > is a positive constant 
related to the volume of A. Since we use this result, we give a version of their proof adapted 
to our situation. 

Lemma 2.1 (Lewis, Pule, and de Smedt [10]). Takee > 0, and take a real L 1 -function 
of positive type v : H u — > R with v(0) > (1-4) ■ There exist a subset A min C H u such that 
for each open box A, with A min C A, the following inequality holds 

£ 'fc-^-^+^fl-)" 1 , (2-4) 

l<i<j<n 

for all n > 2, and each set of n distinct points {xi, . . . , x n } C A. 

Proof. By Bochner's theorem [8] a function of positive type v defines a positive-definite 
quadratic form (•,•)„ on the space of bounded measures on R" : 



(pi, H*)v = v{x-y) d//i(x)d// 2 (y). 



Hence, by the Cauchy-Schwarz inequality this quadratic form satisfies 

, v > \(H2,Vi)v\ 2 

Applying this inequality to the measures 

n 

d^i(x) = ^25 Xt (x)dx and dfx 2 (x) = XB(A,h)(x)dx, 

i=l 

where Xb(a,H)(x) is the indicator of the set B(A,h) = {x E H u : dist(x, A) < h}, and 
Xi e A, i — 1, . . . , n, n > 1, we arrive at the following estimate: 

t v(x '- Xl)> -BS^&y (2 - 5) 

where 

A h A (y) = / dxv(x-y). 

J B(K,h) 

A lower bound for the numerator in the r.h.s. of (2.5) is found using that for all y G A one 
has : 

\A h A (y)-v(0)\< A h A (y)- [ dxv(x-y) + f dx \v(x)\ 

J\x-v\<h J\x\>h 



< 6(h), (2.6) 
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where 6(h) is given by: 



6(h) 



S\x\ 



2 I dx \v(x)\. (2.7) 

'\x\>h 

The last step in (2.6) is valid since for every y e A, there is a ball with the radius h and 
the centre at y lying inside B(A, h). Since v is an ^-function, 6(h) (2.7) converges to zero 
in the limit h — > oo. This yields 

> «(0) - 6(h) > 0, 

for all y G A and for h large enough. The last bound is valid since -0(0) > 0, and since 
6(h) > (2.7) becomes arbitrarily small for h large enough. The numerator in the r.h.s. 
of (2.5) is therefore bounded from below as follows : 



>n 2 (v(0)-6(h)y. 



(2.8) 



An upper bound for the denominator in the r.h.s. of (2.5) is based on the following estimate 



dy A\(y) - Vv(0) 



B(A,h) 



dy A\(y) - [dy A h A (y) 

) J A 



< : 

] B{A,h) 

<vo\(B(A,h)\A)\\v\\ L i + V6(h), 



dy A\(y) - Vv(0) 

(2.9) 



see (2.6). Hence, using the estimates (2.8) and (2.9) in (2.5), one gets 



E 



V(Xi 



Xj) > 



n 



(v(0)-6(h)Y 



V v(0) + ||u|| L ivol(S(A, h)\A)/V + 6(h) ' 



Notice that for any smooth-shaped A the factor vol(i?(A, h)\A)/V is of the order 0(h/L), 
and vanishes in the limit V — > oo for fixed h. Since 



(v(0)-6(h)f 



v(0) + ||v|| L ivol(B(A, h)\A)/V + 6(h) 



> v( Q ) ~ ||u||livo1(S(A, h)\A)/V - 36(h), 



we can choose h large enough, such that 6(h) /v(0) < e/4, and then take A min that 
||u|| L ivol(S(A, h)\A) /Vv(0) < e/4, to obtain the estimate : 



S ^v(x i -x j ) > yv(0)(l-e), 



for all boxes A with A m ,„ C A. 



□ 



As it follows from the proof, this result holds for more general shapes of A, the only 
condition is that vol(5(A, h)\A)/V tends to zero in the limit V — > oo for some fixed h > 0. 
In fact, A tends to H u in the sense of Van Hove [9] is enough. 
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Remark 2.2. From the proof it also follows that the value of e is defined by A min , and 
that increasing the latter we can make e as small as we want. This means that after the 
thermodynamic limit one can put e = 0. 

In [10], it was also shown that the constants A = -0(0) (1 — e) and B = v(0)/2 in 
Lemma 2.1 are optimal for this type of pair-potentials. This is based on the following 
argument: suppose there exists a series of positive constants {Ai}i converging for / — > oo 
to a better superstability constant A > -0(0), i.e. such that for all / 

v(x i -x j )>yA l , (2.10) 

l<i<j<n 1 

for all finite sets of distinct points {x±, . . . ,x n } C A;. Since we can choose e > small 
enough that A — e > v(0), there exists li(e) such that A\ > A — e/2 for all / > l\(e). On 
the other hand, for all A/ large enough, / > /2(e), we get 

/ dxv(x-y) < v(0) + e/2, 

J A t 

uniformly in y G W. Then by integration of both sides of (2.10) over A™, for I > 
max(li(e), 12(e)) we get the estimates : 

nv(0) + n( - n ~ 1) (v(0)+^)>^ [ dxi...dx„ ^v(xi- Xj )>^r(A-€/2). 

Since above the n > 1 is arbitrary, these estimates imply that A — e < "0(0), which 
is in contradiction to the hypothesis, hence yielding the optimality of the constants in 
Lemma 2.1. 



2.2 Thermodynamics of the Mean-Field Bose Gas 

In this section we briefly review the properties of the so-called mean-field Bose gas (some- 
times also called the imperfect Bose gas), an exactly solvable model of Bosons [13-19] (for 
an extended review see [6]), which will play the role of a reference system. The mean-field 
Bose gas is defined by the local Hamiltonians 

Ht = T£ + ±Ni (2.11) 

The kinetic energy operator T A (1.2) is identical to the one of the fully interacting system 
(1.1), but the interaction term (1.3) is replaced by a kind of mean-field interaction term. 
The physical relation of the model (2.11) to our system (1.1) lies in the fact that it is the 
Van der Waals limit of the fully interacting system (1.1) [2]. In that case, the constant A 
has to be chosen equal to the long-range part of the interaction (1.3), i.e. A = gv(0). 



8 



The explicit solution for the gapless case A = of this model can be found at several 
places [6,13-19]. Here, we focus on the case with non- vanishing gap. The grand-canonical 
pressure function at inverse temperature (3 and chemical potential p is defined as 

Trjr s denotes the trace over the boson Fock space Tb- Below, we develop explicit expres- 
sions for the pressure and the particle densities in the thermodynamic limit V — > oo. 

Lemma 2.3 (Thermodynamic Functions). The grand- canonical pressure p A (/3, p) = 
liniAPA a(A p) °f the mean-field Bose gas (2.11) exists for all (5 > 0, p e R and is given 
by the Legendre transform: 

p A ((3, p) = sup (pp-/ A A (/3,p)), (2.12) 

where the canonical free energy f A ((3,p) at inverse temperature (3 and density p is given by 

f£(P,p) = f PA W,p) + ^p 2 , (2.13) 
f p ' A (/3,p) is the free energy of the perfect Bose gas with gap A (1.2). 

Proof. The thermodynamic pressure of the perfect Bose gas is given by: 
p p ' A ((3, p) = lirnp£ A (/3, p) = lim -L Tr^ fl e^"^ 

A A jjv 

which implies that in order to be well defined, p must be bounded from above: p < —A, 



i.e. 



1 oo 

™o=0 {n k }k^o 



The canonical free energy / p,A (/3,p), is the Legendre transform of p p,A ((3, p), defined only 
for p < -A, 

f PA ((3,p) = sup (pp -p pA (/?,p)) . (2.14) 

H<-A 

By direct calculation one finds expression (2.13) for the free energy of the mean- field model 
(2.11) at temperature (3 and density p as 

/ A (/5,p) = lim-^ 7 lnTr <) e-^, 

where Tr_(„) denotes the trace over the Hilbert space of symmetrised functions for 

n = [pV\ (integer part of pV) Bosons. Since on this space the mean-field interaction term 
is constant, we immediately find : 

limf A [H A x }((3,p) = lim f\[T A ](j3, p) + \p\ 
9 



The pressure of the mean-field gas, is again the Legendre transform of f^(f3,p), yielding 
formula (2.12), well defined for all /jgR. □ 

Theorem 2.4 (Pressure of the Mean-Field Bose Gas). The grand-canonical pressure 
of the mean-field Bose Gas (2.11) is explicitly given by 

pt(M J^-^ *" ^-A + V(A-A); „ 15) 

\ (p + Af/2X + p p (P,-A), for ji > -A + Ap p (/3, —A), 

where p p (f3,p), and p p (/3, p) are respectively the pressure, and the total density of the 
perfect Bose gas; p^~°\(3,p) is the pressure of the mean-field Bose gas without gap. 

Proof. The formula (2.15) is found by working out explicitly the Legendre transforms 
in (2.12)-(2.14), and using the properties of the perfect Bose gas. Let p = p\({3,p) be 
solution of the equation: 

for a given p and /3, where the right-hand side is the expectation value of the total density 
in the grand-canonical Gibbs state for the perfect Bose gas model Tfc. Denote the limiting 
solution by 

fi((3,p) = lim# A (/3,p), 
We have p((3, p) < -A, if p < p p ((3, -A), and p((3, p) = -A, if p > p p ((3, -A), hence 

PA ( pp((3,p)-p p ((3,mp)), if P<P P (A-A); 

\-pA-p p (p,-A), if p>p p (f3,-A). 

This is the explicit expression for (2.14). The thermodynamic potentials such as the 
pressure and the particle density of the free Bose gas with gap A are the same as for 
the gapless perfect Bose gas, but only for the values of p, < —A. At p, — —A, there 
is degeneracy of the densities and BEC occurs. Since for A > the critical density 
Pc' A (P) = P P (Pi~A) is finite in all dimensions v > 1, the condensation takes place in 
all dimensions whereas in the gapless case p P '^ A ~°\f3) = p P {(3) < oo only for dimensions 
v > 2. Hence condensation occurs only in three or more dimensions at (3 < oo. 

Recall now the expression for the canonical free energy of the mean-field Bose gas, 
Eq. (2.13). Using the expression for the free energy of the perfect Bose gas f p ' A (/3,p) 
(2.16) derived above, one finds for d p fx((3,p), 

, fA(R . j d p f p (P,p)+Xp, if p<p p (P,-A); 

I -A + Ap, if p > p {(3, -A). 
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By virtue of (2.17) and (2.12), one gets the expression for the the mean-field Bose gas 
pressure: 

Pa (P, p) = SU P {PP ~ f\(P, P)) 

pp(P,p)-ft(P,p(P,p)), for /i<-A + Ap p (A-A); 

M/i + A)/A-/ A A (A(/i + A)/A), for p > -A + Ap p (/3, -A), 

where p((3, p) is the solution of d p f p (0, p(/3, p)) + Ap(/3,p) = p as a function of p < 
-A + Ap p 03, -A) and /3. Since by (2.16) 

/ p ' A (/3, (p + A)/A) = - p P (/3, -A), 

for p + A > Ap p (/3, —A), where p p (/3, —A) is the pressure of the free Bose gas, we use 
the expression (2.13) for the free energy of the mean-field gas to find (2.15), that proves 
Theorem 2.4. □ 

Theorem 2.5. Considering the mean-field Bose gas (2.11), we derive the following ex- 
pressions for the densities in the thermodynamic limit. The total grand- canonical density 
is given by 

^HK^' ^ ' i - 4 + W - i|i (2.18) 
1 (p + A) /A, for /j, > —A + \p p (B, -A). 

TTie zero-mode condensate density is given by 

A f 0, /or p<-A + Ap p (/?,-A); 

P °' 1 (p + A)/A-p P 03,-A), /or p > -A + \p p ((3, -A). 

TTie limit of the expectation value (N^) h a^((3, p)/V 2 is given by 

lim _1 (J v IV (,,,) = ( orW)'. ^ - A + ^ -A); 
A a/h„^,w | ( M + A )2/v, /or ^ > -A + \p p (P, -A). 

ifere p^~°\(3 , p) is the density for the gapless mean-field gas, i.e. for A = in Eq. (2.11). 

Proof. These quantities are derived using that the pressure (2.15) is a convex function 
of respectively A, p and A. By the Griffith lemma [6, Appendix C], the order of the 
thermodynamic limit and the corresponding derivative can be interchanged, which gives : 

Po,\{P>v) = n P t?{ n o)h£ x (P,p) = iipdAPA^Ap) = <9 A p A (/3,p); 
Px(P,p) = li ™y( N *)H£ x (P,p) = lim9 M p A A (/3,p) = <9 M p A (/3,p); 
\im^- 2 (N 2 A ) Hix ((3,p) = lim-2d A p A A (/?,p) = -2d A p A (A p). 
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By virtue of (2.15) of Theorem 2.4 these imply the explicit expressions (2.18)-(2.20) of 
Theorem 2.5. □ 

Taking the limit A j 0, we recover the usual expressions for the mean- field Bose gas 
(2.11) with vanishing gap, in particular the expression for the zero-mode condensate density 
in dimensions v > 2, 

(A=o),„ x / °> for//< Apf(/3); 

Po,x W)H I , P(R . . ^ . Pf „. (2-21) 

3 Proofs of the Main Results 

The main idea of the proof of the condensation for the systems (1.1) is to estimate their 
Bose condensate from below by the condensate of a particularly chosen reference system 
for which one can compute the amount of the condensate explicitly. Thus, a judicious 
choice of this reference system is a subtle point of our proof. 

Since we consider superstable systems, i.e. systems where the grand-canonical pressure 
is defined for all values of the chemical potential, it seems to be natural to choose a reference 
system which is also superstable. This immediately rules out the perfect Bose gas (1.2) as a 
reference system, since its pressure is only well defined for p < —A. Choosing the reference 
systems within the class of mean-field Bose gases (cf. section 2.2), which are indeed well- 
known superstable systems, seems therefore a good choice. The reference systems that we 
consider are mean-field Bose systems which are close enough to the Van der Waals limit 
of the fully interacting system (1.1). Apart from the use of a reference system, the proof 
is based on various convexity properties of the thermodynamic functions. In particular it 
is based on the following lemma. 

Lemma 3.1. The zero-mode condensate density p^ g (/3,p) in the thermodynamic limit 
of grand- canonical Gibbs states of interacting system (1.1) with a superstable two-body 
potentials v satisfying the conditions (a) and (b), has the following lower bound : 

- ^p p (b, -a) - f m- 

Here p g A ~°\(3, p) denotes the total density of the interacting gas without gap (1.1). p P {(3, — A) 
refers to the total density of the perfect Bose gas at the inverse temperature (3 and the chem- 
ical potential p = —A, p P {f3) is the critical density of the perfect Bose gas. The bound is 
valid for values p > gv(0)p P (fl), and dimensions v > 2. 

Proof. The pressure p\[H^] of systems with a gap in the kinetic energy spectrum (1.2) 
and any stable interaction is an increasing convex function of the parameter A > 0. Since 
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by Theorem 2.5 the condensate density Pq 9 (P, //) is the derivative of the corresponding 
pressure with respect to A, the convexity property yields a lower bound for the zero-mode 
density (N /V) h a : 

V {No)H t - A • (3 ' 2) 

Now we use a reference system to get the lower bound on the condensate. This reference 
system is a mean-field Bose gas (2.11), defined by the local Hamiltonian (cf. section 2.2) 

Ht,A = T£ - ^N A + g^N 2 A . (3.3) 

Therefore, we fix the mean-field Bose gas (2.11) interaction parameter by taking A = gA, 
where g > is the coupling constant (cf. (1.1)) and A = t>(0)(l — e) is the optimal 
superstability constant (2.2) associated with the two-body interaction (1.3) of the full 
model (1.1). By virtue of the same convexity property as in (3.2), the difference of the 
pressures between the reference Bose system (3.3) with gap and without gap, is bounded 
from below by the condensate density for the reference mean-field gas without gap, i.e. 



PA[Hj 9 , A ] - Pa[H?J] 1 

A -V (JVo, Cf' (3 ' 4) 



Adding the inequality (3.4) to the lower bound on the condensate density of the full system 
(3.2), we introduce the reference system (3.3) in our estimate: 

^( N o)n tg > ^°><->-X {P^ H LJ -Pa[<= a 0) ] -Pa[H&] + PA[H[ A g =0) ]) . (3.5) 

Hence, the condensate density of the interacting model (1.1) with gap (A > 0) is bounded 
from below by the condensate density of the mean-field model (3.3) without gap (A = 0), 
and a correction term proportional to 1/A containing the pressure differences between the 
full system and the reference system. 

These pressure differences will be estimated using the Bogoliubov convexity inequality 
[6, Appendix D]. Applied to the grand-canonical pressures of the mean-field reference Bose 
gas (3.3) and the full model (1.1), it gives 

f (W£) Hta < p A [<,J - p A [<] < 9 - (W£) Htg A , (3.6) 

for any A > 0. Here the operator W A is the difference between the interactions of the fully 
interacting and the mean-field Bose gases : W A = U\ — AN\/2V . Then by virtue of (3.5) 
and (3.6) we get : 
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Now our task is to estimate the two expectation values of in (3.7). An upper 
bound on (W^/V) h a in (3.7) can be found using the properties of the pair-potential v 

and the Gibbs states of the reference system (3.3). Expressed in terms of the creation and 
annihilation operators on A*, we get for (W A /V) h a ^: 

qeA* fceA* fc'eA* 



f^ 4 E E ( a i a k>a{a k ) H A g A . 



2V 2 

fceA* fc'eA* 



Exploiting the mode by mode gauge invariance of the Gibbs states of the mean-field Bose 
gas (3.3), and rewriting the above expression in terms of the occupation-number operators 
Nk = a\ak we arrive at 



v«)^ A = E E (*(°) + *(* - k ') - A ) ( N * N «)»t, A 

k V 

^(o)E(w 2 W„,. 4 + w^„J- 



(3.8) 



fc 



Since by condition (b) : v(0) > v(k) > 0, the coefficients in the first sum of the r.h.s. of 
(3.8) are bounded as 

\ (v(0) + v(k - k') — A) < v(0) - A/2. 

From the second sum in the r.h.s. of (3.8), we retain only the quadratic zero-mode term, 
by the Cauchy-Schwarz inequality we have, 

This yields the following upper bound for {W^/V) h a : 

V {W ^ )h L,a ^ 2V 2 Wh£ 9iA - zvt( N °)H tgtA (3-9) 

The expectation values appearing in the r.h.s. of (3.9) can be calculated exactly, ap- 
plying Theorem 2.5. They give in the thermodynamic limit the upper bound : (v(0) — 

A/2)pfJf3^) 2 -v(0)ptAP^) 2 / 2 - 

The other unknown term in (3.7) is (Wy^)„(A=o). It can be estimated using the su- 

A, 9 

perstability (2.3) of the interaction (1.3) by the tuning the interaction parameter of 
the mean-field reference Bose gas (3.3) to be equal to the constant A in the superstability 
criterion (2.3), which gives the estimate from below : 

^«>„(a=o) > -^(N A ) H(A=0) . (3.10) 

V A, g V JJ A, 9 
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This, in particular, justifies our choice of the parameter A = gA specifying the reference 
system (3.3). Using now (3.9) and (3.10) in (3.7) one finds in the thermodynamic limit the 
following lower bound for the condensate density p^ g (/3,p) 

- f (B Pg A =°H(3,p) + (v(0) - A/2)p* A {(3,pf) . 

The lower bound (3.1) now follows from the explicit expressions (Theorem 2.5) for the 
total density and the condensate density of the mean-field Bose gas with gap, and from the 
well-known expression for the condensate density in the gapless mean-field model (2.21) 
for p > gApc(/3), i.e. in the regime where p^ g ~A > 0. 

In the last step to (3.1) we use the optimal superstability constants for continuous L l - 
potentials of positive type (cf. Lemma 2.1), i.e. we put A = (1 — e)-fi(O), and B = v(0)/2. 
This gives the expression for the lower bound in the form (3.1), since by Remark 2.2 we 
can put e = after the thermodynamic limit. □ 

Notice that the lower bound (3.1) contains the term p g A ~°\(3,p), i.e. the total density 
of the fully interacting gas without gap. It is not explicitly known as a function of (3 and p. 
However it is always finite, and it can be viewed as a reference parameter. Using a slightly 
modified reference system, an alternative lower bound can be derived which consists only 
of explicitly known functions related to the perfect Bose gas. 

Lemma 3.2. The zero-mode condensate density p A g ((3,p) in the thermodynamic limit of 
the grand- canonical Gibbs states of interacting systems (1.1) with superstable two-body 
potential v satisfying conditions (a) and (b), has the following alternative lower bound: 

a (a \ ^ ty + 9v(0) gv(0) P(a A , 2 P( „s 
P ^ A ~ 2gv(0) + ^A P (A " A) " Pc iP) 

2p + 2A + gmm + mpP{ ^ 



2Av(0) 

p p (P,—A) refers to the total density of the perfect Bose gas at inverse temperature (3 and 
chemical potential p = —A, and Pc(P) is the critical density of the perfect Bose gas. The 
bound is valid for all values p > gv(0)p^((3), and dimensions v > 2. 

Proof. The proof is completely analogous to the proof of Lemma 3.1. But now we use 
the alternative reference system : 

Hi 9 , A ,B = T£ -pN^ + g [^Nl - BN^ , (3.13) 

which compared to the first reference system (3.3), contains an extra interaction term. 
Since the term —gBN\ is linear in the total number operator, it corresponds to a shift in the 
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chemical potential. Again, the constants A and B coincide with the optimal superstability 
values (2.2) for the pair-potential v of the full system (1.1), where g > 0. 

First, we derive a bound similar to the one of (3.7) in the proof of Lemma 3.1. Now 
one gets : 

- " i {v {w *' B) »t,^ ~ • <*•"> 

where W A ' B = U A - AN\/2V + BN A . The expectation values in the r.h.s. of (3.14) can 
be estimated analogously to (3.9) and (3.10). This yields for the upper bound : 

1 IW A < B \ <T 2 ^(°) ~ A /Ar2\ ^(°)/AM2 B IN\ 

> h L,a, b ^ 2V* {Na)h L^b - 2V^ {N ° )H t,A, B + v {Nk)h -^b- ( 3 - 15 ) 



For the lower bound of (W A ' B ) h (a=o) we use again the superstability of the interaction U 



'A, 9 



A 



(1.3), and the fact that according to the superstability criterion (2.3) we can take for A 
and B their optimal values. This gives: 

hw A > B ) H( ^ > 0. (3.16) 

V A, 9 

The explicit formula (3.12) now follows if one introduces (3.15) and (3.16) into (3.14), 
using the explicit expressions for the densities of the mean-field Bose gas (section 2.2), and 
for the optimal values of A and B (section 2.1), and finally taking taking e = after the 
thermodynamic limit, see Remark 2.2. □ 

It should be remarked that one can hardly compare the bound given in Lemma 3.1 
with the one in Lemma 3.2, and hence to express an opinion which of them yields the best 
result. However, as the latter bound is known explicitly, it can be used to make numerical 
estimates of the condensate density and of the minimal gap as functions of the various 
parameters involved. For this we refer to ref. [11] and section 4. 

Instead, we proceed now with the proof of our Theorem 1.1, based on the lower bound 
derived in Lemma 3.1. We prove that the condensate density of the full model (1.1) is 
strictly positive in the domain jj, > gv(0)p P (j3) if the gap is large enough. 

Proof of Theorem 1.1. Consider the bound from Lemma 3.1, 



+ eMw, -Af - ^l P r o) (A /*) - ^V(A -A). 



(3.17) 



Fix the inverse temperature (3 and take the chemical potential /i compatible with the 
condition of the theorem. This ensures the positivity of the first term in the r.h.s. of 
(3.11)). Now take \i such that 
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for some arbitrarily chosen r] > 0. This yields a lower bound for the first line in the r.h.s. 
of (3.17). 

The expression on the second line can be absolutely bounded by rj using A large enough 
and the fact that liniA->oo P P (f3, —A) = 0. This gives : 



,f( V(A^-#»r»»/')- !L+ij ' 



2A 



2A 



-A) 



< V, 



for all A larger than some minimal gap: A > A min , which exists for the fixed \i > 
gv(0)Pc(P)- 

Collecting these two estimates, we obtain that for a fixed temperature and rj > one 
can find \i and A large enough such that : 



d£ g (P,n)>ri>o, 



proving the condensation. 



□ 



Similarly one can prove the existence of the zero-mode condensation on the basis of the 
bound found in Lemma 3.2. 



4 Discussion 

So far, we are concentrated on the case of dimensions v > 2, however, the result of 
Theorem 1.1 can be extended to dimensions v — 1 or v — 2. A lower bound for the 
condensate density p^ g (P, p) for v < 2 as in Lemma 3.1 or Lemma 3.2 is derived in a similar 
way. It requires slightly modified convexity arguments (3.2)-(3.4). Since the free Bose gas 
in dimensions v < 2 shows only condensation in the case of non-vanishing gap (1.2), one 
has to consider in (3.2)-(3.4) the pressure differences in the form p\[H^] — p\[H^ }, for 
some A > 0, with < A < A, instead of A = 0. This yields the substitution in (3.11) 
of Pq^~a \fi \ p) an d Pg A ~°\[3,ij) by p Ao gA ((3,p) and pf°([3,fi). The bounds derived in this 
way are valid for all dimensions v > 1, and lead to similar conclusions as in Theorem 1.1. 
Hence, in one and two dimensional interacting Bose gases with large enough gap (1.2), the 
zero-mode Bose-Einstein condensation is also proved. Notice that this is in contrast to the 
Bogoliubov-Hohenberg theorem [12, 20] which yields the absence of BEC for translation 
invariant continuous Bose systems without gap for dimensions v < 2. 

We use the lower bound (3.12), which can be computed explicitly as a function of the 
different parameters, to visualise our estimates. In Figure 1, the dependence of the lower 
bound on the temperature is indicated. The lines shown on the (/z, A)-graph indicate 
domains where the lower bound (3.12) is positive, i.e. above each of these curves we have 
BEC. The dashed curve is the threshold for condensation calculated for inverse temperature 
(3 = 0.05, and the solid line is the threshold at inverse temperature (3 = 0.1. Clearly, for 
higher values of f3, the condensation occurs for smaller gaps and for smaller values of /i, 
i.e. at lower densities. 
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Figure 1: (fi, A)-graph with temperature dependence. 

To get an idea of the phase diagram of our model, on Figure 2 we present a family 
of thresholds as a function of the gap value A, the plotted curves are the thresholds for 
A = 0.6, 0.8, and 2. The dotted line is the line fi = gv(0)p^(f3), it indicates the border of 
validity of our estimates. As above, this family is calculated by equalising the lower bound 
(3.12) to zero. Notice that to get the real phase diagram one has to do this for the left 
hand side of (3.1) and not for the lower bound. 

Considering the high density (large /i) regime, the lower bound (3.12) can be written 

as 

" f 1 1 w°Kmirv,-x)))+ow, 



0(0) Vfl A V 2 

which means that in order to have a positive lower bound for p^ g ((3,/i), we need 

A>g(^-+v(0)p P (P,-^)], (4-1) 



for high values of fi. This means that for non-zero interaction (1.3) (g > 0), there is a 
non-zero lower bound on the gap width. Notice that the minorant (4.1) is proportional to 
the coupling constant g > 0. 

If we now choose for the two-body interaction a family of Van der Waals scaled pair 
potentials, i.e. we substitute 

v{x-y)*-+\ v v{\{x-y)), (4.2) 
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Figure 2: (//, (3) phase-diagram with A-dependence. 

A > 0, in the expression for the interaction term (1.3), then we find that for a fixed A > 0, 
the condition (4.1) is satisfied if A is chosen small enough and for low enough temperatures. 
This is easily seen as follows: substitution (4.2) amounts to substituting A"t>(0) for v(Q) 
and leaving v(0) in (4.1) unchanged. Therefore the r.h.s. of (4.1) can be made smaller than 
any A > 0, by choosing A small enough for large /3 such that p p {[3, — A) gets small. Hence, 
we recover the result of Buffet, de Smedt and Pule [2] about the stability of Bose-Einstein 
condensation in the weakly interacting Bose gases with the Van der Waals scaled potentials 
and a non-zero one-particle spectral gap. 

If v(0) tends to infinity, then the condition (4.1) can not be satisfied for any finite 
gap. In this case our estimate becomes a triviality. This behaviour is compatible with the 
observation that for the hard core continuous Bose gas in a scaled external field there is no 
condensation with a macroscopic occupation of any level of the one-particle Hamiltonian 
[21]. 

Finally we remark that our results are for continuous homogeneous systems. The only 
assumptions we make are the gap in the one-particle excitations spectrum (1.2) and the 
superstability conditions on the pair-potential v (conditions (a) and (&)). Various other 
interesting exact results on Bose condensation are known, e.g. for Bose systems with a 
family of Van der Waals potentials [2,22], for models with truncated interactions [6,23,24], 
or for Bose lattice models with hard core interaction and at half-filling [25]. Only recently, a 
proof of BEC is found for the trapped interacting gases [1], i.e. for inhomogeneous systems, 
in the so-called Gross-Pitaevskii limit for particle interactions. 
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